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ABSTRACT 

The performance of an integrated GPS/inertial attitude determination system is 
investigated using a linear covariance analysis. The principles of GPS interferometry 
are reviewed, and the major error sources of both interferometers and gyroscopes are 
discussed and modeled. A new figure of merit, Attitude Dilution of Precision (ADOP), 
is defined for two possible GPS attitude determination methods, namely single differ- 
ence and double difference interferometry. Based on this figure of merit, a satellite 
selection scheme is proposed. The performance of the integrated GPS/inertial attitude 
determination system is determined using a linear covariance analysis. Based on this 
analysis, it is concluded that the baseline errors (i.e., knowledge of the GPS interfer- 
ometer baseline relative to the vehicle coordinate system) are the limiting factor in 
system performance. By reducing baseline errors, it should be possible to use lower 
quality gyroscopes without significantly reducing performance. For the cases consid- 
ered, single difference interferometry is only marginally better than double difference 
interferometry. Finally, the performance of the system is found to be relatively insensi- 
tive to the satellite selection technique. 
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Chapter 1 


Introduction 

1.1 Motivation 

Attitude determination is a vital function for all spacecraft. Accuracy requirements 
are often very strict to accommodate precise pointing applications. The sensors cur- 
rently used for attitude determination include star trackers, sun sensors, and gyro- 
scopes. A common configuration is comprised of high accuracy gyroscopes and star 
trackers. The gyroscopes are used to provide continuous data, but because all gyro- 
scopes are subject to drift, star trackers or other sensors are needed as an inertial refer- 
ence to realign the gyroscopes periodically. GPS interferometry may be used as that 
inertial attitude reference. With the increasing use of GPS for navigation, interferome- 
try may be added with only a small incremental cost and change, as opposed to adding 
completely separate system, such as a star tracker. 

GPS interferometry uses carrier phase measurements to determine vehicle attitude. 
Phase measurements are made at two or more antennas, and then differenced to pro- 
duce a precise measurement of the direction to the GPS satellite relative to the baseline 
connecting the antennas. This measurement is then combined with the vehicle and 
GPS satellite positions to produce an inertial attitude measurement. 

The use of GPS interferometry for space vehicle attitude determination is not a 
new idea. Tests of airborne and ground based GPS interferometers have indicated that 
accuracy sufficient for many space applications is possible with the current technology 
[9, 18]. Trimble and Loral recently announced that they will team to develop a space- 
qualified receiver capable of performing vehicle attitude determination. Also, GPS is 
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being incorporated into the attitude determination system for the International Space 
Station. 

Stand-alone GPS attitude determination systems have been developed for marine 
and other earth-based applications [5, 9], Typically, attitude measurements are only 
available at about 1 Hz, and no measurement is available if the GPS receiver is unable 
to maintain track on the signals of a sufficient number of satellites. This is not accept- 
able for space vehicle navigation, where continuous attitude estimates are often 
required. If the GPS interferometer is coupled with a gyroscopic attitude reference sys- 
tem, then the latter provides continuous data output and stability during short periods 
of high dynamics. Integration of GPS interferometers with gyroscopes has been pro- 
posed for a variety of applications, including autonomous space vehicle navigation 
[24], and vehicle attitude determination [2, 9], including spacecraft [6, 23]. Because 
GPS provides an inertial reference at every measurement, it is not necessary for the 
gyroscopes to maintain accuracy for long periods of time. Therefore, a GPS/inertial 
system could make use of smaller, cheaper gyroscopes than are currently used in most 
space applications. 

An overview of the principles of GPS interferometry is presented in the next sec- 
tion, followed by a discussion of the major performance factors which are addressed in 
this thesis. Finally, the goals and organization of this thesis are described. 

1.2 GPS Interferometry 

As shown in Figure 1.1, attitude determination using GPS is accomplished by 
measuring the difference in carrier phase as the signal arrives at two separate antennas. 
The total phase difference consists of a fractional part, A<|>, which is often referred to as 
the phase difference, and a number of whole wavelengths, n, called the integer ambi- 
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Figure 1.1 GPS Signal Arrives at Two Antennas 

guity. Only the fractional part is measured because the integer ambiguity is not observ- 
able. Once computed, the integer ambiguity acts as a bias. It will be discussed further 
in Chapter 2. The total phase difference is related to the range difference from the sat- 
ellite to each of the antennas (leg A of the triangle shown) by 

A<|> , 

A = (n + 2 ^)*. (>•») 

where X. is the known carrier wavelength. The phase difference measurement is called 
a single difference (SD). 

The orientation of the antenna baseline in inertial space relative to the line of sight 
(LOS) to a satellite is the angle 0, which is related to the phase difference and ambigu- 
ity by 

A«b X 

cosG = (n + — ) - = u b - u R (1.2) 

where is a unit vector along the baseline, and Ur is a unit vector along the LOS to 
the satellite. In general, the resolution of 0 improves as the baseline gets longer. 
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Determination of the three dimensional baseline orientation in space is equivalent 
to the calculation of the coordinates of the baseline vector, b. Since there are three 
unknowns, three independent SD measurements are required to completely determine 
the orientation. The phase difference measurements are not sensitive to vehicle roll 
about the axis of the baseline; therefore, at least two baselines are needed to determine 
the three dimensional attitude of a vehicle in space. 

Interferometry may also be accomplished using double differences (DD) or triple 
differences (TD). The DD technique involves differencing phase difference measure- 
ments from two satellites (see Figure 1.2). Notice that the notation for a single differ- 
ence is A, which has one vertex on top (for the one satellite in the sky), and two 
vertices on the bottom (for the two antennas). The notation for a double difference, V, 
has two vertices on top, representing the two satellites used for the measurement. With 
the DD technique, some common mode receiver errors cancel; however, an additional 
satellite is required, since four satellites are needed to get three independent measure- 
ments. Triple differencing involves differencing DD measurements from one time 
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period to the next. TD measurements can be used to eliminate the integer ambiguity 
[14, 16]. 

1.3 Performance Factors 

One important aspect of implementing GPS interferometry is satellite selection 
[5]. From low earth orbit, there are often as many as 10 or 1 1 satellites in view. All-in- 
view-processing, which is becoming more common as advances in receivers are made, 
requires much more additional hardware for attitude determination than for position- 
ing, and is not currently affordable for most applications. 

The optimal satellite geometry for attitude determination is not the same as that for 
positioning [13]. Therefore, the dilution of precision (DOP) metrics defined for posi- 
tioning, such as geometric dilution of precision (GDOP) and position dilution of preci- 
sion (PDOP), are not valid for attitude determination. Although satellite selection 
schemes for pointing applications have been presented [11, 14], no instantaneous fig- 
ure of merit for three dimensional attitude determination has been defined. 

Another important factor is the measurement type. Double differences seem to be 
the most common type [3, 5, 9, 22], possibly because they are traditionally used in sur- 
veying and other differential GPS applications. In these applications, long baselines 
are employed, and DD are used in order to eliminate the offset between the receiver 
clocks which exists because separate receivers are used to make measurements at each 
antenna. In attitude determination, all of the phase measurements are made in refer- 
ence to a common oscillator in a single receiver, so there is no receiver clock error [ 1 , 
18]. Double differences can be useful in attitude determination, because there may be 
differences in the electrical path length between each antenna and the receiver [1,5, 
18]. These path delays are common to all satellites being tracked, and therefore cancel 
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out in double difference measurements. However, there are disadvantages in using 
double differences, because the measurement noise level is increased by combining 
the measurements, and more satellites are required to obtain the same number of mea- 
surements. 

1.4 Thesis Objectives 

This thesis will develop a model and linear covariance simulation of a coupled 
GPS/inertial system for a space vehicle in low earth orbit. Performance characteristics 
of the coupled system will be investigated for both good and poor quality gyroscopes. 
The effects of using single difference instead of double difference measurements will 
also be investigated. A definition of attitude dilution of precision (ADOP) will be 
developed, and the performance of various satellite selection schemes will be com- 
pared. 

The thesis is organized as follows. Chapter 2 discusses the major error sources for 
both GPS interferometers and gyroscopes. Models for the errors included in the simu- 
lation are presented and explained. Chapter 3 presents details of the satellite selection 
routines. ADOP is derived, and sample calculations are presented. In Chapter 4, the 
details of a linear covariance analysis are provided. The implementation used in this 
thesis is described, and a flow chart of the simulation is presented. Performance results 
are presented and discussed in Chapter 5. In Chapter 6, a summary of the thesis is pre- 
sented, and conclusions are drawn from the results. Suggestions for future study are 
also discussed. 
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Chapter 2 


Error Sources and Modeling 

Both GPS interferometers and gyroscopes are subject to errors from a number of 
sources. These errors can affect the performance of the GPS/inertial attitude determi- 
nation system. For this reason, it is important to understand the sources of the errors. 
This chapter first describes the gyroscope/interferometry configuration with respect to 
the vehicle and the flight attitude. Subsequently, it discusses the principal errors that 
affect GPS/inertial systems, and lays out the technique that was used to model them 
for this analysis. 

2.1 Gyroscope/Interferometer Configuration 

The mounting configuration between the interferometer and the gyroscopes was 

modeled as a rigid attachment, disregarding any flexures that may exist between the 
two. The baselines of the interferometer are orthogonally aligned along the x- and y- 
axes of the vehicle body. 

The LVLH coordinate frame is defined in a non-traditional manner for this discus- 
sion with the z-axis aligned along the inertial position vector, as illustrated in Figure 
2.1 

2.2 GPS Error Sources 

GPS interferometer errors are discussed in detail below, along with their effect on 
attitude measurement. Any compensation techniques that can be used to reduce the 
errors are also described. 
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Figure 2.1 LVLH Coordinate Frame Definition 

2.2.1 Phase Error 

Phase error is the error made by the receiver when it computes the phase angle of 
the incoming signal. Phase error depends mostly on the quality of the receiver and fre- 
quency of the output, but may also be affected by high vehicle dynamics if the vehicle 
dynamics are not within the tracking loop bandwidth. If the GPS tracking loops are 
inertially aided (as in a coupled system), then normal spacecraft dynamics will not 
contribute to phase error. 

2.2.2 Multipath 

Multipath is caused by reflecting surfaces near the antenna which cause the signal 
to arrive at the antenna via more than one path, or worse yet, via only a non-direct path 
(see Figure 2.2). Multipath errors tend to be less severe in attitude determination appli- 
cations, because the measurement is a difference between two closely spaced anten- 


20 



nas. Therefore, the majority of multipath signals are “common mode,” i.e., they affect 
both antennas equally. “Differential mode” multipath does not cancel, and tends to 
increase as the baselines get longer. 



Figure 2.2 Signal Multipath 


Multipath errors can be limited in several ways, including careful antenna place- 
ment on the vehicle, appropriate choice of antenna gain pattern and axial ratio, the use 
of choke rings, and coupling the GPS with a gyroscopic attitude reference system. 
Antennas should not be mounted near obvious sources of multipath. Most multipath 
enters the antenna at low elevation angles, so an antenna with a gain pattern which has 
low gain in this area will be less susceptible to multipath. Choke rings can also be used 
to mask signals which would enter the antenna at low elevation angles. 

Axial ratio refers to an antenna’s sensitivity to signal polarization. All broadcast 
GPS signals have the same polarization, which becomes reversed if the signal reflects 
off another surface. By attenuating signals with incorrect polarization, the effect of 
most multipath signals will be minimized. Of course, signals which have reflected off 
of two surfaces would not be rejected in this way, but will usually be significantly 
attenuated by the multiple reflections. 
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When GPS signals corrupted with multipath are used in the attitude solution, it 
appears that an antenna or baseline has moved or rotated from its true orientation. If 
the multipath signal is being reflected off of the surface on which the antennas are 
mounted, then the apparent rotation will vary slowly as a sinusoid which is related to 
the repeating geometry of the satellite line of sight with respect to the antenna, and the 
reflecting object. This type of error is very difficult to distinguish from actual vehicle 
motion. Fortunately, most multipath of this type is common mode multipath, and can- 
cels out in the phase difference measurement. The remaining multipath can possibly 
be estimated and calibrated, as previously described [4], 

An object which passes too close to the antenna array, for example, the Shuttle 
Remote Manipulator System, or another spacecraft performing docking or rendezvous 
procedures, could cause either common mode or differential mode multipath. This 
type of multipath may begin and end abruptly as the object moves past the antenna 
array, causing severe errors in the attitude solution. However, if the GPS interferome- 
ter is coupled with gyroscopes, then the latter will provide a record of vehicle rota- 
tions, and it is possible to discount signals which indicate a sudden vehicle rotation 
which is inconsistent with the gyroscope state. 

2.2.3 Flexures 

There are two types of flexures which can limit the accuracy of a GPS attitude 
determination system: flexures of the antenna baseline, and structural flexures between 
the baseline and the navigation base of the vehicle. The first case is illustrated in Fig- 
ure 2.3. It is possible to estimate the change in baseline length [2], but the change in 
orientation is difficult to distinguish from vehicle rotation. This type of flexure can be 
minimized by keeping the baseline as short as possible. 
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Original Baseline 

Flexed Baseline 
Measured 

movement Baseline 

Figure 2.3 Flexures Cause Apparent Change in Baseline Length and Orientation 

Flexures between the antenna baseline and the navigation base cause changes in 
the orientation of the baseline with respect to the orientation of the vehicle. Both types 
of flexures can be caused by thermal expansion which may occur when the spacecraft 
passes in and out of the earth’s shadow, and by any stresses which occur during 
launch. 

2.2.4 Propagation Errors 

GPS signals are refracted and delayed as they pass through the atmosphere, result- 
ing in propagation errors. The two parts of the atmosphere which cause propagation 
errors are the troposphere, which extends up to about 100 km, and the ionosphere, 
which extends from 100 km to 1000 km above the earth’s surface. Spacecraft in low 
earth orbit, which are considered in this thesis, operate in the ionosphere. Thus all GPS 
signals received by such spacecraft are subject to ionospheric refraction and retarda- 
tion. In addition, measurements through the earth’s limb are also subject to tropo- 
spheric errors. However, in attitude determination, the signal travels very nearly 
identical paths to reach each antenna, and virtually all atmospheric delays cancel out 
when the phase difference is computed. Signal refraction does not cancel, but the total 
bending through the atmosphere is only about 0.003° [25], and this effect is negligible 
compared to other error sources. 
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2.2.5 Clock Errors 

Clock errors exist in both the GPS space segment and user clocks. When the phase 
difference between signals from the same satellite is calculated, the bulk of these 
errors cancel out. The small satellite clock error which remains due to the different 
transmission times of the signals which arrive at the separate antennas is on the order 
of 10 seconds for the highly accurate satellite clocks, and is not significant com- 
pared to the other error sources. Since all phase measurements are made against a 
common oscillator, user clock errors also cancel out. 

2.2.6 Path Delay 

Temperature differences in the cabling and electronics result in a difference in 
electrical path lengths between each antenna and the receiver. The path delay is some- 
times called the “apparent user clock bias” because it resembles the bias seen in Dif- 
ferential GPS applications when the phase measurements are made on different 
receivers (and therefore not relative to a common oscillator.) The path delay may be 
eliminated by using the double differencing (DD) technique. 

2.2.7 Integer Ambiguity and Cycle Slip 

The integer ambiguity represents the number of whole wavelengths in the phase 

difference at two antennas. It is not directly measurable. There is one integer ambigu- 
ity for each satellite per baseline. In other words, if three satellites are being tracked 
using three antennas (and two baselines) there will be six integer ambiguities to be 
determined. Once an integer ambiguity has been determined, it is maintained as long 
as the tracking loop maintains lock on the satellite signal; a cycle slip occurs if the 
tracking loop loses lock on the GPS signal. If undetected, a cycle slip causes the previ- 
ously calculated integer ambiguity to be invalid. An incorrect integer ambiguity causes 
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a bias-type error in the attitude measurement. Because of its quantized nature, the inte- 
ger ambiguity cannot be estimated in a linear Kalman filter. Search algorithms have 
been developed which exhaustively check the possibilities with varying degrees of 
efficiency. These algorithms usually require extensive computation, the use of extra 
satellites, and/or special receiver characteristics [7, 10]. At least one commercial man- 
ufacturer uses baselines shorter than the LI wavelength (X, = 19 cm) which eliminates 
the ambiguity [5]; but accuracy of this configuration is limited by the extremely short 
baseline. For high accuracy applications, a better solution is to couple the interferome- 
ter with a gyroscopic attitude reference system. When an attitude estimate is provided 
by gyroscopes, the integer ambiguity may be instantaneously determined, or highly 
constrained, depending on the accuracy of the estimate. 



Figure 2.4 Angular Effect of the Integer Ambiguity 

For example, consider a simple case where n=0 or n=l (illustrated in Figure 2.4). 
5/2 is the angular resolution required to resolve the integer ambiguity. It depends on 
the baseline length, b, and the direction to the satellite, 0. For all values of 0, 5 and b 
are related by 





For baselines of 1 meter, an attitude estimate accurate to ±5° is sufficient to uniquely 
determine the integer ambiguity. Then the integer ambiguity may be calculated at each 
measurement, rather than being a stored value. In this case, cycle slips would not 
occur. Alternately, since inertial systems are very reliable for short periods during high 
dynamics, they may be used both to detect the initial cycle slip, and to maintain an atti- 
tude measurement until GPS signal lock can be regained. 

2.2.8 Satellite Geometry 

Poor satellite geometry can degrade the accuracy of the GPS attitude measure- 
ment just as it degrades position and velocity measurements. However, the satellite 
geometry requirements are slightly different for attitude determination. As in the posi- 
tioning problem, it is desirable that the satellites be spread out on the sky; but in addi- 
tion the line of sight (LOS) to each satellite should be as perpendicular as possible to 
each baseline. When the satellites are not spread out, there tends to be a bias in the atti- 
tude solution [5, 12] (just as there would be for the position solution). As stated in 
Equation 1 .2, the value of cos0 is calculated for the solution. Because the slope of cos0 
is steepest around 0 = 90°, the best angular resolution of 0 is obtained in that region. 
If double differences are used, the difference of the LOS vectors should be perpendic- 
ular to the baselines. The most precise measurements can be obtained if the satellite 
selection routine is applied separately for each baseline, rather than for the array as a 
whole [5], In this case, different satellites may be selected for each baseline, increasing 
the number of receiver channels required and the complexity of the processing soft- 
ware. 

At this point, no instantaneous figure of merit for attitude determination has been 
reported in the literature. Later in this thesis, a definition of ADOP will be presented, 
along with further discussion of the merits of various satellite selection routines. 
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2.3 Gyroscope Errors 

The following describes the major error sources for ring laser gyros (RLGs), the 
type under consideration in this thesis. 

2.3.1 Bias 

All gyros are subject to bias drift rate, which causes a drift in the attitude solution 
over time. When the gyroscopes are coupled with a GPS interferometer, GPS measure- 
ments are available about once a second. Therefore, as long as unacceptable drift does 
not occur between measurements, a very large bias can be tolerated. In addition, nearly 
all of the bias is typically a random constant, which can easily be estimated in a Kal- 
man filter if precise measurements are available. 

2.3.2 Scale Factor Error 

The angular output of a gyroscope varies proportionally to the angular input. How- 
ever, the constant of proportionality is not known precisely, causing the slope of the 
output vs. input line to differ from the design slope, as shown in Figure 2.5. The differ- 
ence between the true slope and the design slope is the scale factor error. 



2.3.3 Misalignment 

The gyro misalignments are also called “non-orthogonalities.” Because of manu- 
facturing limitations, the three single-degree-of-freedom gyroscopes required for an 
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attitude reference system cannot be mounted perfectly orthogonal, resulting in the mis- 
alignments. There are two misalignments for each gyro; for example, the x-axis gyro 
has one misalignment about the y-axis, and another about the z-axis. Therefore, there 
are a total of six misalignments. The misalignments may also be estimated, provided 
that the vehicle performs maneuvers which make the misalignments observable, and 
GPS or other reference measurements are available. 

2.3.4 Angle Random Walk 

Noise on the angular rate output causes the attitude solution to wander or “random 
walk.” The filtering process can estimate part of the resulting attitude error, but cannot, 
of course, predict how the error will random walk. Therefore, angle random walk will 
directly limit the length of time that gyroscope accuracy can be maintained without 
GPS measurements or other alignments. 

2.4 GPS Error Models 

Statistical error models are used to capture the effects of the errors in the linear 
covariance analysis. The phase error, including receiver noise, is modeled as zero 
mean white noise. The baseline flexures and path delays are both modeled as zero 
mean band-limited noise. Physically, that means that the baselines cannot flex with 
infinite frequency. It also means that the flexures are temporally correlated. The same 
is true of the electrical path length changes which cause the path delays. 

This type of random process is adequately approximated by an exponentially cor- 
related random variable (ECRV) [19], and has the autocorrelation function 

_|T| 

<P XX (t) = a 2 e T (2.2) 
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where x is the time difference variable, and T is the time constant. In the remainder of 
this thesis, x is used to denote the time constant of an ECRV. The time constant for 
both the baseline flexures and the path delays is tied to the heating and cooling which 
is caused by the vehicle moving in and out of the sunlight. 

Multipath errors are also modeled as ECRVs [17, 20]; the multipath time constant 
is the time it takes for the LOS to the satellite to “significantly” change direction, or 
move through about 30° [25], The multipath error on a measurement depends on the 
direction to the satellite; since the satellites are selected to be spread out in the sky, the 
multipath error from each satellite is modeled as independent. The multipath error 
from one satellite to each of the orthogonal baselines may be independent or highly 
correlated, depending on the local reflective environment. Thus the total multipath on 
a particular measurement, q t , is modeled as 

( 2 - 3 ) 

where is the common multipath error, and ^ is the multipath error which is indepen- 
dent, or different between the two baselines. 

Propagation errors and satellite clock errors are considered negligible (as dis- 
cussed in the previous chapter), and ignored. It is assumed that a batch processing 
solution resolves the initial integer ambiguity (at system start-up) [10]. Once an atti- 
tude estimate exists, integer ambiguity can be calculated at each measurement using 
the gyroscopes, so that cycle slip is prevented. The interferometer/navigation base 
flexures were not modeled, but they may not be negligible. Any flexure of this type 
directly contributes to the attitude errors. 
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2.5 Gyroscope Error Models 

Gyroscope bias, scale factor, and misalignments are all modeled as random con- 
stants. Although long term instability does exist in each error source, it is not signifi- 
cant for the durations considered in this thesis. 

Because the angular rate is integrated to produce attitude, the noise on the angular 
rate is also integrated. The effect of the noise is modeled as a random walk, which is 
valid to the extent that the noise is actually uncorrelated (white). 
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Chapter 3 


Satellite Selection 

3.1 Introduction 

Poor satellite selection can contribute to the attitude measurement errors. Chapter 
2 described in qualitative terms the recommended satellite geometry. In order to 
describe the satellite geometry quantitatively, it is necessary to define a numerical fig- 
ure of merit. The literature discusses two possible ways to do this: AZDOP [13], and 
Pitch and Yaw Sensitivities (abbreviated PYS) [14], AZDOP is a one dimensional fig- 
ure of merit, i.e., it applies when only the heading is being measured. PYS is a two 
dimensional routine. While they provide useful information, both of these methods are 
inadequate for the three dimensional attitude problem presented here. This chapter 
presents the derivation of a three dimensional figure of merit, attitude dilution of preci- 
sion (ADOP), for satellite selection. Sample calculations and geometries are 
described. An expanded, three dimensional version of Pitch and Yaw Sensitivities is 
also presented. GDOP was used in the simulation for comparison purposes, and is 
described in Appendix A. 

3.2 Attitude Dilution of Precision 

The effect of satellite geometry on measurement error is often described using 
dilution of precision factors. For example, if each measurement in a set has a standard 
deviation given by CT 0 , then the standard deviation of a parameter of interest, 0, will be 
related to 0 O by [16] 

a = DOPo 0 (3.1) 

DOP factors provide convenient, instantaneous figures of merit for selecting satellites. 
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The derivation of the DOP factor for three dimensional attitude (ADOP) used in 
this analysis is similar to the derivation of GDOP described in Leick [16]. First, 
Leick’s definition of the design matrix is applied to the attitude determination problem. 
This yields the design matrix 


Parameters 


A = 


Roll 

5 ^n° bS| 


5sn° bs > 


Pitch 

v ft obs. 
9pitch 


Yaw 


(etc.) 


0 

ST 

1 

to 


The observables may either be single differences (SD) or double differences (DD). 
(The sensitivities which appear in the design matrix are derived in the following sec- 
tions.) Then 


(A t A)" = 


a «t>0 a 4>v 
*e<t> °e 


( 3 . 2 ) 


where (|>, 0 , and \|/ represent roll, pitch, and yaw. The off-diagonal terms represent the 
cross correlations, so, for example, represents p (j)e <J (|) (T e , where p^ e is the correla- 
tion coefficient between roll and pitch errors. Assuming that the measurement noises 
are uncorrelated, ADOP is related to the design matrix by 


ADOP = Jtr ( A t A) 


-i 


( 3 . 3 ) 
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which means that 


ADOP = fif+cl + cl (3.4) 

If the measurement noises are correlated, then ADOP is related to the design matrix by 

ADOP = 7tr(A T R"'A) _l (3.5) 

The positive definite weighting matrix, R, which is related to the correlation of the 
measurement errors, has the form 


R = 


2 1 1 
1 2 1 
1 1 2 


(3.6) 


given that the double differences are all related back to a single satellite. These double 
differences can be computed in alternate manners; but the corresponding form of R 
changes such that ADOP does not depend on the double difference definition. 


3.2.1 Variations and Extensions of ADOP 

Just as PDOP is extended to GDOP by including time, it may be desirable to 

include other parameters in ADOP. In particular, in the SD case the path delays may be 
included. Baseline coordinates could also be included. The number of other parame- 
ters that are included is limited by the number of independent measurements available, 
and hence the number of satellites used, since A*A must have full rank for an inverse 
to exist. 


3.2.2 Single Difference Sensitivities 

All sensitivities were derived from the single difference measurement equation, 
which (omitting measurement noise sources) is 

SD = (u^^bf + p, (3.7) 


33 



In Equation (3.7) the superscripts I and B indicate a vector coordinatized in the earth 
centered inertial frame, and the vehicle body frame, respectively, and the subscript i 
indicates which baseline was used for the measurement. The path delay is represented 
by p, and Cg is the current body to inertial transformation matrix. 

The first step in the derivation is to express Cg, the vehicle attitude, in terms of the 
parameters of interest: roll, (J), pitch, 0, and yaw, \| f. Roll, pitch, and yaw are defined to 
be ordered, fixed-axis rotations. Therefore, Cg is related to (J), 0, and \j/ by 



c\|i — s\|/ 0 

c0 0 s0 

10 o’ 

R z (\|/)R y (0)R x (<t>) = 

svp c\p 0 
.0 0 1 

0 1 0 
— s0 0 c0 

0 c<|> — s<() 
0 s<j) c<{) 


where c\j/ means cos\| /, s0 means sin0 and so forth. When Cg is expanded 


C 


i 

B 


cy • C 0 C\)/ • S 0 • s<|> - S\|/ • c<|> c\|/ ■ S 0 • C<(> + S\]/ ■ S<t> 

SV|/ C 0 sy • S 0 • S<t> + C\J/ C<|) SX4/ ■ S 0 • C<(> — C\|/ - S<() 

— S 0 C 0 • S<t> C 0 • C<J) 


The baselines used in this thesis were (bf ) T = [1 0 o| and (b®) 1 = [0 1 6] 
Expanding the single difference for the first baseline gives 


SD = [ U x U y U z] C E 


1 

0 

0 


+ P. 


= U x • cy • C 0 + u y • sv • C 0 - u z • S 0 + p, (3 g) 

Then the sensitivities for the measurements taken across baseline 1 can be determined 
by taking partial derivatives of Equation (3.8), so that 
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li SD 

5§ SD 

^-SD 

dy 

_d_ 

dp, 

j_ 
dp 2 


= o 


= -u x • CA|/ S0 - U y Sl|/ • S0 - u 2 C0 
= — U x • C0 • SVj/ + U y • C0 • 


C\|/ 


SD = 1 


SD = 0 


Likewise, the expansion of the SD for the second baseline yields 


( 3 . 9 ) 


SD = [u x u y U Z ]C[ 


+ P 2 


= — u x s\|/ • c<|> + u y • C\|/ - c<> + u x • C\|/ • S0 • S<t> 
+ U y • SVj/ ■ S0 • S<|) + u z • C0 • S<t> + p 2 


( 3 . 10 ) 


and the corresponding sensitivities are 


J^SD = u x • s\|f • s<(> — u y c\|/ • s(|> + u x • c\|/ ■ s0 • c<|) + u y svj/ s0 • c<|> 

+ u z • c0 • c<}> 

^SD = u cv|/ s<f> • c0 + u • sv)/ ■ s<{> • c0 - u z ■ s<J» • s0 
ou y 

^SD = -u x • c<j) • c\y - u y • c<t> • sip - u x • s0 • s<J) sip + u y ■ s0 • s<|) • c\|/ 


dp 


dp 


-SD = 0 

i 

-SD = 1 


( 3 . 11 ) 


3.2.3 Double Difference Sensitivities 

In the DD case, the measurement equation (omitting measurement noise sources) 


is 
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DD = (Ur, — uj^ 2 ) T C^b® (3.12) 

Thus, the sensitivities for the DD case can easily be obtained by substituting the differ- 


ence of the LOS vectors, (u RI - u R2 ) , for the LOS to a single satellite, u R , in the sin- 
gle difference equations. Therefore, the sensitivities for baseline 1 are 

l-DD = 0 

d<p 

JqDD = -(u xl -u x2 ) CM/ s© - ( Uy , — u y2 ) S\|/ S0- (u zl -u z2 ) c0 

^-DD = -(u xl -u x2 ) c0 s\j/+ (u yl -u y2 ) • c0 ■ c\|/ 
y (3.13) 

and the sensitivities for baseline 2 are 

^DD = (u xl — u x2 ) sy s<|> — (u yl -u y2 ) • c\j/ • s<j> + (u xl -u x2 )c\|/s0c<|> 

+ (u y , - u y2 ) ■ svj/ • s0 • c<}> + (u zl - u z2 ) c0 c<> 

^DD = (u xl -u x2 ) c\|/ • s<|> • c0 + (u y ,-u y2 ) s\|/ - s<|> c0 

“( U zl “ U Z2) S<}) S© 

^DD = - (u xl - u x2 ) c<t> ■ cm/ - (u yl -u y2 )c<t>- s\|f- (u xl - u x2 ) s0 • s<|> • s\|/ 

+ (Uy, - u y2 ) S0S<1>C\1/ (314) 

Since the path delays are not observable in the DD case, they are not included in 
ADOP. 

3.2.4 Sample Calculations of ADOP 

This section describes the calculation of ADOP for the satellite geometry shown in 

Figure 3.1, when single difference measurements are used. Subsequently, the position 
of satellite 3 is varied to show how ADOP changes with the changes in satellite geom- 
etry. The path delays are not included in these calculations. 
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For the sample calculations, it is assumed that the body coordinate frame is aligned 



2 

Figure 3.1 Satellite Geometry 

with the inertial frame, so that 

<j) = e = y = 0 (3.15) 

The unit LOS vectors for the satellite geometry shown are 

Uri = [0 0 l] 

y £ 2 = [i o o] 

UR3 = [o 1 o] 

Using Equations (3.9) and (3.11) to calculate the sensitivities. 
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The first two rows of A are the sensitivities of the single difference using satellite 1 
and baselines 1 and 2, respectively. Equation (3.16) shows that measurements taken 
from the satellite overhead are sensitive to either roll or pitch, depending on which 
baseline is used. Measurements from either of the satellites on the horizon are sensi- 
tive to yaw only. Then 


(A t A)"' 


1.0 0.0 0.0 
0.0 1.0 0.0 
0.0 0.0 0.5 


and for this satellite geometry 


ADOP = J~( A t A)"‘ = 1.58 

Because the baselines are orthogonal, ADOP does not vary when the azimuth of satel- 
lite 3 changes. Figure 3.2 shows how ADOP varies with the elevation of satellite 3. 



Figure 3.2 ADOP Changes With Satellite Geometry 
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3.3 Pitch and Yaw Sensitivities 

This algorithm is based on a two dimensional algorithm presented previously [14], 
That algorithm included equations for the elevation and azimuth sensitivities of a 
baseline which was nearly aligned with a local level frame, such that the antennas 
were always pointed away from the earth. Those equations have been adapted for 
space applications, where the orientation of the antenna array is completely arbitrary, 
and renamed the pitch and yaw sensitivities. The derivation of the equations used in 
this thesis from those presented in the reference is shown in Appendix B. 



Figure 3.3 Pitch and Yaw Sensitivities 


In the context of this algorithm, “pitch” and “yaw” refer to baseline movement 
about a nominal position, rather than spacecraft attitude. In Figure 3.3, pitch is rotation 
out of the antenna array plane, or about the c-axis, while yaw is rotation in the antenna 
array plane, about the a-axis. 

“Pitch sensitivity” is the sensitivity of a particular SD or DD to baseline orienta- 
tion in the pitch direction. “Yaw sensitivity” is the sensitivity to baseline orientation in 
the yaw direction. Both sensitivities are defined in terms of the spherical coordinates 
of the satellite with respect to the baseline in question: 
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PS = sirup 


( 3 . 17 ) 


YS = cosrpsinG ( 3 . 18 ) 

If DD measurements are used, then the spherical coordinates of the difference of the 
LOS vectors, rather than the LOS to an individual satellite, are used when computing 
the sensitivities. 

The pitch or yaw sensitivity of a group of four satellites (or three double differ- 
ences) is calculated by taking the sum-squared of the individual sensitivities. Thus, 
each group has two total sensitivities, rather than a single figure of merit. Groups are 
rated against one another by comparing the minimum (or worst) sensitivities; the 
group with the largest (best) minimum sensitivity is ranked the highest. 
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Chapter 4 


Linear Covariance Analysis 

A linear covariance analysis computes the error statistics of a system. Therefore, it 
can be used as a tool used to predict system performance. This chapter begins with an 
overview of linear covariance analysis. Then, the simulation used for this thesis is 
described, and a flow chart is presented. Finally, the equations used to implement the 
linear covariance analysis in this simulation are described. 

4.1 Overview 

The linear covariance analysis is based on assumptions that the dynamics of the 
system are linear or can be linearized, that an estimate of the system state is main- 
tained, that the errors in that estimate can be statistically described, and that these sec- 
ond-order statistics can be propagated in time. 

The dynamics of the system state are described by 

x = Fx + w (4.1) 

where x is the state vector, F is the dynamics matrix, and w is Gaussian white noise. 
The analysis in this thesis focuses on the system state at discrete points in time, t k , 
t k+ |, ...; in that case the discrete dynamics can be described by a difference equation 
[ 8 ], 

X(t k+I ) = 4>(t k+1 ,t k )x(t k ) +w k (4.2) 

where <J>(t k+1 , l k) is the state transition matrix for the time step from t k to t k+1 , and w k 
is the discretized noise. w k is related to w by 
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(4.3) 


W k = J 0(t k+1 ,T)w(T)dT 

•k 

While the actual value of the system state is not known, an estimate, denoted by x, 
is maintained. The errors in the estimate, called the state errors, are defined as 

5x k = X k -X k (4.4) 

These errors are described by their second-order statistics which together are called 
the covariance matrix, E. The covariance of the state errors is defined as 

E k = (5x k 5 x k > (4.5) 

where ( ) denotes the expectation operation. 

If the initial value of E is known, then the knowledge of the state dynamics and any 
measurements that are performed can be used to calculate future values of E, or in 
other words, to provide a covariance analysis. The covariance analysis has two basic 
parts: (1) the process propagation, and (2) the measurements. The process is how the 
errors change due to the state dynamics alone. When the covariance matrix is known at 
time t^, then it can be propagated to a future time t^+i using [8] 

E(t k+I ) = 0(t k+ „t k )E(t k )cD T (t k+1 ,t k ) +Q k (4.6) 

where Q k is the discrete process noise matrix, and is also called the noise covariance 
matrix. Q k is related to the process noise spectral density matrix, Q, by 

<k + i 

Qk = J <I>(t k+l ,T)Q(T)d> T (t k+1 ,T)dX (4.7) 

‘k 

and to the discretized noise by 

Qk = <™k ' w k > (4.8) 
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The process propagation evolves the covariance matrix to the times at which the mea- 
surements occur. 

The second major part of the linear covariance analysis is the measurement incor- 
poration. Measurements are taken in order to reduce the state errors. It is assumed that 
each group of simultaneous measurements, 2 , can be modeled as a linear combination 
of the state plus some Gaussian noise. The relationship between the measurements and 
the state is given by 

Z k = H k x k + v k (4.9) 

The observation matrix, H k , and the measurement noise, v k , as well the method used 
to incorporate the new measurement with the previous estimate, determine the effect 
of a particular measurement on the error covariance matrix. Although there are many 
approaches to incorporating measurements, the one which is most common is the Kal- 
man filter equation. If the state estimate is updated in this manner, then the new covari- 
ance matrix is related to the old by [8] 

E + = E“-E~H k (H k E~H k + R k ) 'H k E _ (4.10) 

where the measurement noise matrix, R, is defined as: 

R k = <y k -v^) (4.ii) 

E _ and E + are used to designate the covariance matrix before and after the measure- 
ment has been incorporated, respectively. 

If the measurement errors, v k , are not correlated with each other, then R k will have 
the diagonal form 


43 



(4.12) 



< V I V I> < V I V 2> < V ! V 3> ••• 


a* 0 0 ... 

R. = 

<v 2 v,) <v 2 v 2 > <v 2 v 3 > ... 


0 a 2 0 

k 

<v 3 V|) <v 3 v 2 > <v 3 v 3 > ... 


0 0 a] 


L ’ ' ’ * * ’ 


... 


Then the measurements can be incorporated sequentially, rather than in a batch. In this 
case, each individual measurement is modeled as 

z k = b k * k + v k (4.13) 

where h^ is the sensitivity vector, v k is the measurement noise, and (v k ) = a^. Each 
measurement can then be incorporated using 

E + = E _ -E'b k (b k E _ b k + a^) _l h k E" ( 4 . 14 ) 

Sequential incorporation is sometimes preferred, because the measurements do not 
have to be taken at exactly the same time, and it is possible to incorporate some of the 
measurements, even if others are not available. There is also a computational advan- 
tage, in that there is no need to take a matrix inverse in Equation (4. 14). 

4.2 Simulation 

The linear covariance analysis developed for this thesis includes a simulation of a 
spacecraft in low earth orbit. Orbital conditions typical for a shuttle or station type 
spacecraft were chosen. Attitude estimates are provided by a gyroscopic attitude refer- 
ence system, which are updated with GPS carrier phase difference measurements. A 
flow chart for the simulation used for this thesis is shown in Figure 4.1. 

Also included is a simulation of the GPS constellation, which approximates the 
satellites’ orbits as circular, with an altitude of approximately 26,600 km. The mod- 
elled constellation contains 24 satellites, arranged in 6 rings of 4 satellites each. The 
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Figure 4.1 Simulation Flow Chart 
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rings are equally spaced in argument of latitude, and at an inclination of 55° . Table 4. 1 
gives the details of the position of each satellite in its ring. 


Table 4.1 GPS Satellite Constellation 


Ring 

Number 

Longitude of 
the Ascending 
Node (deg) 

True Anomaly for Each Satellite in the 
Ring (deg) 

1 

0.0 

190.96 

220.48 

330.17 

83.58 

2 

60.0 

249.90 

352.12 

25.25 

124.10 

3 

120.0 

286.20 

48.94 



4 

180.0 

312.30 

340.93 

87.06 

209.81 

5 

240.0 

11.90 

1 10.76 

143.88 

246.11 

6 

300.0 

52.42 

165.83 

275.52 

305.04 


4.3 Implementation 

In this section, the equations used to implement the linear covariance analysis are 
presented. Included are the equations for the error models discussed in Chapter 2, and 
the details of the covariance propagation. The complete error state vector for the SD 
case is x T = [ £ T $ T y T 5b| SbJ £ T p} p] ]. Each element of the error state is 
described in Table 4. 1. A more detailed description of each error was given in Chapter 
2. 

In the DD case, the path delays are not included in the state. The following sec- 
tions describe the error state dynamics and the measurement incorporation. 

4.3.1 Process Dynamics 

The state transition matrix is used to propagate the error state covariance matrix 
according to Equation (4.6). Since there is no closed form solution for the state transi- 
tion matrix, it is calculated using the approximation [8] 
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Table 4.2 Elements of The Error State Vector 


Symbol 

Description 

Number of 
Elements 

HL 

Attitude Error 

3 

£ 

Gyro Bias Error 

3 


Gyro Scale Factor Error 

3 

* 

Gyro Misalignment 

6 

6b, 

Baseline 1 Error 

3 

8b 2 

Baseline 2 Error 

3 

l 

Multipath Error 

12 

Pi 

Path Delay 1 

1 

P2 

Path Delay 2 

1 


0 = I + FAt + -2j-F 2 + ... (H.O.T.) (4.15) 


which is truncated at the second-order term. It is assumed that F is, or can be approxi- 
mated as, constant over At. Sufficient accuracy is ensured by comparing the fourth 
term in the sequence with the 0 approximated as the sum of the first three. A warning 
is printed if the ratio of any element in the fourth term to its corresponding element in 
the approximated 0 exceeds one one-thousandth. 


4.3. 1.1 Attitude Errors and Gyroscope Dynamics 

The error in the attitude solution that is produced by the gyroscopes is determined 

by the various gyro errors: bias, e, scale factor, p, the misalignments, y, and the white 
noise T), which produces angle random walk. The dynamics of the gyro process, coor- 
dinatized in the body frame are 


47 




vj/ B = (— cpg x ij/ B ) + e B + 


e 1 

X 

O 

o 

1 


o CD y 0 

J3 B + 

0 0 CD Z 



-(0, co 0 0 0 0 

z y 

0 0 0) z -0) x 0 0 

0 0 0 0 -CD CD 

y 


7® + T) B (4.16) 


where \j/® is the attitude error vector, CDg is the body rotation rate with respect to the 
inertial frame, and cd x , co y , and co z are the components of CDg. The first four terms of 
Equation (4.16) appear directly in the dynamics matrix. Since T| is a noise term, it 
appears in the spectral density matrix, Q, rather than in the dynamics matrix. The com- 
plete dynamics and noise matrices are shown in Appendix C. 


4.3. 1.2 Baseline Errors, Path Delays, and Multipath 

Since each baseline coordinate error, each path delay, and the multipath error asso- 
ciated with each satellite is modeled as an ECRV, the dynamics of these errors are 
independent and have the same form. The dynamics for a general ECRV, x, are [8] 

1 

X = - - X + V (4.17) 

where X is the time constant, and v is white noise. The discretized covariance of v over 
a time step At depends on the variance of the random process, o 2 x , the time constant, x, 
and the time step. At. It is given by [8] 


2At > 


1 - e T 


(4.18) 


The complete dynamics matrix and process noise matrix are shown in Appendix C. 


4.3.2 Measurement Updates 

All measurements were incorporated sequentially, using Equation (4.14). As stated 
previously, two types of measurements are used, single differences (SD) and double 
differences (DD). Only one measurement type is used during a particular run. 
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4.3.2. 1 Single Difference Measurements 

When SDs are used, the measurement noises are uncorrelated, and the measure- 
ments can be incorporated directly using the sequential form. The relationship 
between the sensitivity vector and the measurement is shown in Equation (4.13). The 
sensitivity vector can be calculated from Equation (3.7); the complete derivation is 
shown in Appendix D. The nonzero parts of the sensitivity vector are h^ (attitude 
error), h b) or h b2 (baseline errors) depending on whether this measurement was taken 
using baseline 1 or baseline 2, h^ (multipath), and hpj| or h pt j 2 (path delays), again 
depending on the baseline in question. 

The measurement noise term, a 2 , used to incorporate the measurement is simply 
(v 2 >. 

Care must be taken when updating the multipath states. Since each multipath error 
belongs to a certain satellite, when a change in the satellites used for measurement 
occurs, the variance for a new satellite’s multipath must be reset to the initial value, 
and the appropriate cross-correlations must be set to zero. 

4.3.2.2 Double Difference Measurements 

When DD are used, the measurement noises are correlated, and the measurements 

cannot be incorporated directly using the sequential form. However, the measurements 
can still be incorporated sequentially if they are first decorrelated. This is done by 
transforming the measurements so that the new measurement noise matrix has a diag- 
onal form. The details of this process are shown in Appendix D. The measurement 
noise term, or, for each measurement is taken from the diagonal of the new measure- 
ment noise matrix. 

The non-zero parts of the sensitivity vector (before the decorrelating transforma- 
tion) for the DD case are h^ (attitude error), h b) or h b2 (baseline errors) depending on 
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whether this measurement was taken using baseline 1 or baseline 2, and h^ (multi- 
path). As in the SD case, care must be used in regards to the multipath state. 
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Chapter 5 


Results 

The linear covariance simulation presented in the last chapter was used to analyze 
the performance of a GPS/inertial attitude determination system. This chapter summa- 
rizes the results of the simulation runs. A description of the factors which affect the 
system performance is presented, and a nominal run is defined. Finally, the results for 
various performance factors are presented and discussed. 

5.1 Performance Factors and Nominal Conditions 

The performance of a GPS/inertial attitude determination system is affected by 

many factors. In order to understand the effects of each performance factor, they must 
be varied individually. A nominal run was defined which served as the baseline for 
these variations. The factors which were varied in this thesis, along with the range of 
values which were considered, are discussed in the following section. Subsequently, 
the nominal run conditions are described, including those factors which were not var- 
ied in this thesis. 

5.1.1 Performance Factors 

The performance factors which were varied, along with the range of values consid- 
ered, are listed in Table 5.1. The phase error listed is the error on a single phase mea- 
surement. A phase error of 3° is equivalent to a 1.6 mm error in the range to the 
satellite. The multipath correlation refers to the correlation between the multipath from 
one satellite to each of the two baselines. The good quality gyro was selected to have 
characteristics similar to those of the ring laser gyroscope system proposed for the 
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Space Station; the poor quality gyro has significantly larger errors, such as a microme- 
chanical gyroscope might have. 


Table 5.1 Performance Factors 


Performance Factor 

Range of Values Considered 

Measurement Type 

SD, DD 

Satellite Selection 

SD: GDOP, PYS, ADOP, Extended ADOP 
DD: GDOP, PYS, ADOP 

Baseline Flexure Error 

a = 0, 1,5 10 mm 

Multipath Error Level 

0 = 0, 1.9, 3.8 mm 

Multipath Correlation 

0, 50, 100% 

Orbital Inclination 

0°, 51.6°, 90° 

Attitude Hold 

LVLH, Inertial, Maneuvers 

Gyroscope Quality 

Good 

Poor 

Drift Bias (°/hr) 

0.01 

10.0 

Scale Factor Error (ppm) 

2.3 

30.0 

Misalignment (arcsec) 

8.3 

100.0 

Angle Random Walk °! Jhr 

0.0033 

0.5 


5.1.2 Nominal Run Conditions 

Table 5.2 shows the conditions which were defined as the nominal run. The orbital 
altitude and eccentricity were chosen to be typical of a space station type vehicle. The 
same GPS constellation was used for all the runs, so that while the satellites visible to 
the spacecraft may have changed due to the orientation of the spacecraft, the GPS sat- 
ellites were always in the same place at the same time in every run. In general, the 
selection of the GPS measurement rate depends on vehicle dynamics and performance 
requirements, as well as the capability of available equipment. For vehicles with low 
dynamics and long mission periods, such as the Space Station, which is considered in 
this thesis, the steady state accuracy is more pertinent than the short term response. 
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Table 5.2 Nominal Run Conditions 


Factor 

Condition or Value 

Orbital Altitude 

110 nm (204 km) 

Orbital Eccentricity 

0.0 

GPS Constellation 

As Described in Section 4.2 

GPS Measurement Rate 

1 Measurement Every 20 
seconds 

Baseline Configuration 

2 Orthogonal Baselines: 

1 Along the Body X-Axis, 
1 Along the Body Y-Axis 

Baseline Length 

1 m 

Phase Error 

3° 

Path Delays 

a = 2 mm, x = 2657 s 

Baseline Flexure Error 

a = 5 mm, x = 2657 s 

Multipath Error 

a = 1.9 mm, x ~ 770 s 

Multipath Correlation 

50% 

Measurement Type 

SD 

Satellite Selection 

Extended ADOP 

Orbital Inclination 

51.6° 

Attitude Hold 

LVLH 

Gyroscope Quality 

Good 


The effect of the measurement rate on the performance was investigated by conducting 
several runs which differed from the nominal case only by the GPS measurement rate. 
The results of these runs are shown in Figure 5.1. 

It can be seen from the figure that attitude error is estimated more quickly when a 
faster measurement rate is used. However, after about 5000 seconds, the attitude accu- 
racy is nearly the same, regardless of the measurement rate. Based on these results, the 
frequency of GPS measurements was chosen to be only one every 20 seconds. 
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Figure 5.1 Effect of GPS Measurement Rate on Attitude Errors 


although most currently available equipment is capable of producing measurements at 
about 1 Hz. 

The time constants for both the baseline errors and the path delays are equivalent 
to one half of an orbital period. The time constant for each multipath ECRV depends 
on the rate of change of the LOS to the satellite; it computed dynamically in the simu- 
lation. 
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5.2 Results for the Nominal Run 

This section describes the results for the nominal run in detail. Plots of the square 

root of every diagonal element of the state vector error covariance matrix are pre- 
sented and discussed. 

5.2.1 Attitude Error 

Plots of each component of the attitude error, as well as the root sum square (RSS) 
of the components, are shown in Figure 5.2. All of the components of attitude error are 
substantially reduced by a single GPS measurement (see Table 5.3); however, the sys- 
tem has a long settling time, approximately 50,000 s, as indicated on the figure. 

After the first measurement, the z-component of attitude error is smaller than the x- 
and y-components; this difference is due to the baseline configuration. Rotation about 
the body z-axis can be measured using both baselines, while rotation about the body x- 
axis can only be measured using the baseline 2 (which is along the body y-axis), and 
rotation about the body y-axis can only be measured using baseline 1 (which is along 
the body x-axis). However, after about one orbit, the y-component of attitude error is 
much larger than the x- and z-components. This result is related to the nominal attitude 
hold orientation, in which the spacecraft is rotating about its body y-axis at the orbital 


55 



rate. This rotation apparently makes it more difficult to estimate the y-component of 
attitude error. This result will be further investigated in Section 5.3.6. 


Table 5.3 Attitude Errors for the Nominal Case 


Component 

Attitude Error (mrad) 

t = 0 
seconds 

t = 20 
seconds 

t = 6000 
seconds 

t = 86400 
seconds 

X 

17.45 

5.94 

1.20 

0.36 

Y 

17.45 

5.97 

3.45 

1.89 

Z 

17.45 

4.43 

1.20 

0.36 

RSS 

30.23 

9.52 

3.84 

1.96 


5.2.2 Baseline Errors 

The errors in each component of each of the baselines are shown in Figure 5.3. The 
noisy appearance of the baseline errors is due to the fact that they were modeled as 
ECRV’s. Other error sources create a noise floor, below which it is not possible to 
reduce any of the baseline errors. 

The lengths of the baselines (the x-component of baseline 1, and the y-component 
of baseline 2) are estimated more quickly than the other components, because they are 
more easily distinguished from vehicle rotation. The z-component of baseline 1 is 
much more difficult to estimate than any of the other components. Since baseline 1 is 
along the x-axis, a z-component error is equivalent to a rotation about the y-axis, and it 
is difficult for the system to distinguish this baseline flexure from the vehicle rotation, 
which is also about the y-axis. 
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Attitude Error Magnitude (mrad) Z-Axis RMS Error (mrad) Y-Axis RMS Error (mrad) X-Axis RMS Error (mrad) 
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Figure 5.2 Nominal Case Attitude Errors 
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Z-Axis RMS Error (mrad) Y-Axis RMS Error (mrad) X-Axis RMS Error (mrad) 


Baseline 1 RMS Errors 


Baseline 2 RMS Errors 





Figure 5.3 Baseline Errors 
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5.2.3 Multipath Errors 

The total multipath error on a particular measurement was modeled as the combi- 
nation of the multipath which is common to both baselines, and the multipath which is 
independent between the baselines. Thus, the variance of the total multipath error, o t 2 , 
is related to the variances of these components by 

o; = a 2 c + of + 2p cj a c a, (5.1) 

where p ci is the correlation factor between the two multipath components. 

The standard deviation for one of the multipath states for baseline 1 is shown in 
Figure 5.4. The other seven multipath states have very similar settling times and 

2 1 1 1 1 1 i 1 i 1 I 


1.8h 
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0.4- 

0 . 2 - : : - 

qI I I 1 1 1 1 1 1 

01 2345678 

Time (s) x10 4 

Figure 5.4 Nominal Case Multipath Error 

steady state values. As with the baseline errors, the noisy appearance is a result of 
using an ECRV to model this error. There is a noise floor which prevents reduction of 
the multipath error below about 1 .5 mm. 
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An expansion of the first 6000 seconds of Figure 5.4 is shown in Figure 5.5. The 



sharp peaks which are indicated on the graph show when changes were made in the 
satellite selection. These peaks do not extend up to the initial error level because a 
measurement is incorporated immediately after the satellite change, and the post-mea- 
surement value is saved and plotted. The difference between the initial error and the 
level of the peaks indicates the amount of estimation that is possible with a single mea- 
surement. The frequent switching of satellites makes the estimation of the multipath 
error more difficult. 

5.2.4 Path Delays 

The path delays were only estimated when single difference measurements were 
used, since the delays are not observable using double difference measurements. Fig- 
ure 5.6 shows the RMS magnitude of the path delay for baseline 1. As with the other 
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Figure 5.6 Path Delay for Baseline 1 

errors which were modeled as ECRV’s, the path delay is noisy in appearance. There is 
also a noise floor which limits the estimation of the path delay below about 1 mm. 

5.2.5 Gyroscope Errors 

None of the gyroscope errors are estimated very well, because they contribute only 
a small amount of error in the time between each measurement. Figure 5.7a, b, and c 
show the gyro bias, scale factor errors, and misalignments, respectively. Not every 
component of each error is observable in the nominal case, because the vehicle is only 
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Figure 5.7a Gyro Bias Errors 
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Figure 5.7b Gyro Scale Factor Errors 



Time (s) 

Figure 5.7c Gyro Misalignments 
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rotating about the body y-axis. This can be seen in Table 5.4, which shows the values 
of each component at the beginning and end of each run. 


Table 5.4 Gyroscope Errors 


Error 

Component 

a, 

t=0 

seconds 

a, 

t=86400 

seconds 

Gyro Bias 
(deg/hr) 

X 

0.0100 

0.0100 

Y 

0.0100 

0.0069 

Z 

0.0100 

0.0100 

Scale Factor 
Error 
(ppm) 

X 

2.3333 

2.3333 

Y 

2.3333 

2.3313 

Z 

2.3333 

2.3333 

Misalignment 

(arcsec) 

Yxy 

8.3333 

8.3333 

Yxz 

8.3333 

8.3240 

Yyx 

8.3333 

8.3333 

V 

8.3333 

8.3333 

Yzx 

8.3333 

8.3287 

Yzy 

8.3333 

8.3333 


5.3 Results 

In this section, the results of the individual runs for the different performance fac- 
tors are presented. For each case, one factor at a time was varied from the nominal run 
presented in Table 5.2. Eight different performance factors were examined; these were 
summarized in Table 5. 1 . 


63 




5.3.1 Satellite Selection Routine 

The attitude error magnitude was used to compare the performance of the satellite 
selection routines. The results are summarized in Table 5.5. These results show that for 


Table 5.5 Satellite Selection Routine Performance 


Run Description 

Attitude Error 
Magnitude at 
t=6000 s (mrad) 

SD GDOP 

3.84 

SD PYS 

3.85 

SD ADOP 

3.93 

SD Extended ADOP 
(nominal case) 

3.84 

DD GDOP 

3.86 

DD PYS 

3.88 

DD ADOP 

3.87 


applications with long mission durations, the satellite selection is not a very important 
factor in either the single difference or the double difference case. 

5.3.2 Measurement Type 

The factor varied in this case was the measurement type; the system performance 
when double difference measurements were used was compared to the nominal case, 
in which single differences were used. Since the RMS measurement noise on each sin- 
gle difference measurement is lower, by the square root of 2, than the RMS measure- 
ment noise on each double difference measurement, it was thought that if the path 
delays could be estimated well enough to counteract the increased noise, then better 
performance would be obtained using single difference measurements. Another 
advantage of single differences is that they are computationally simpler; double differ- 
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ences are not independent, and must be decorrelated before being incorporated. Also, 
one fewer satellite is needed for unique attitude determination. A disadvantage is that 
two more filter states are needed to estimate the path delays. Table 5.6 shows that there 
is very little difference in system performance due to the measurement type. 


Table 5.6 Measurement Types 


Measurement 

Type 

Attitude Error RMS 
Magnitude at t=6000 sec 
(mrad) 

SD 

3.844 

DD 

3.847 


5.3.3 Baseline Errors 

The performance level is directly affected by how well the baselines are known, as 
shown in Figure 5.8. Large or moderate baseline errors are also responsible for the 
extremely long settling time of the system, as shown in Table 5.7, and Figure 5.8. 


Table 5.7 System Settling Times 


Baseline Errors 
(mm) 

Settling Time 
(seconds) 

0 

10,000 

1 

15,000 

5 

50,000 

10 

80,000+ 


Estimation of the baseline length down to the noise floor is possible even with large 
initial error levels, as shown in Figure 5.9 and Figure 5.10. In each run considered in 
the section, the spacecraft was in the nominal attitude orientation, with the body x-axis 
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Figure 5.8 Attitude Errors due to Different Baseline Error Levels 


pointed in the direction of velocity, and the vehicle rotating about the body y-axis. 
Therefore, baseline 1, which is along the body x-axis, sweeps out 360° during each 
orbit. The body z-axis errors of this baseline are difficult to distinguish from the orbital 
rotation, resulting in poor estimation of this quantity. 


5.3.4 Multipath 

Two facets of the multipath error were varied: the level of error, and the amount of 
correlation between the multipath errors received by the two baselines. The level of 
error was varied by changing the magnitude of the a associated with each multipath 
error state. The amount of correlation was varied by changing the amount of multipath 
which was common to both baselines relative to the amount which was different for 
the two baselines. 
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Z-Axis RMS Error (mrad) Y-Axis RMS Error (mrad) X-Axis RMS Error (mrad) 





Figure 5.9 Baseline 1 Errors 
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Figure 5.10 Baseline 2 Errors 
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5.3.4. 1 Multipath Level 

Multipath level was varied from 0 to 3.8 mm (nominal is 1.9 mm). As expected, 
the larger the multipath errors are, the less accurately the attitude can be estimated. 
(See Table 5.8.) However, there is only a small improvement in attitude error between 
the zero and nominal multipath cases. This indicates that the nominal multipath level 
modeled here does not contribute significantly to the attitude errors. 


Table 5.8 Multipath Level 


Multipath Level 
a, mm 

Attitude Errors 
at t=6000 seconds 
(mrad) 

X 

D 

Z 

RSS 

0 

1.14 

3.40 

1.15 

3.76 

1.9 

1.20 

3.45 

1.20 

3.84 

3.8 

1.30 

3.56 

1.30 

4.01 


5.3.4.2 Multipath Correlation 

The amount of multipath correlation between the baselines was varied from 0- 
100%, using the relations 


af = ( 1 - p) a 


(5.2) 


and 


°l = P<J 2 (5.3) 

where is the nominal variance of the multipath error, and are the variances 
of the independent and common multipath errors, respectively, and p is the correlation. 

Slightly better attitude estimates are obtained when the multipath is highly corre- 
lated, as shown in Table 5.9. In this case, there are only four multipath processes 
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which are being estimated (one for each satellite), instead of eight or twelve as in the 
non-correlated or partially correlated cases. However, the table also shows that there is 
no significant difference in the accuracy of the attitude estimate due to the amount of 
multipath correlation. 


Table 5.9 Multipath Correlation 


Multipath 

Correlation 

Attitude Errors 
at t=6000 seconds 
(mrad) 

X 

Y 

Z 

RSS 

0% 

1.20 

3.45 

1.20 

3.85 

50% 

1.20 

3.45 

1.20 

3.84 

100% 

1.19 

3.45 

1.19 

3.83 


5.3.5 Inclination 

In a highly inclined orbit, it was thought that there may be fewer satellites visible 
than in an equatorial orbit or an orbit with the nominal inclination of 51.6°. If fewer 
satellites are available, then the performance of a GPS-based system could be limited 


in high inclination orbits. However, as Table 5.10 shows 

Table 5.10 Orbital Inclination 


Orbit Inclination 

Attitude Error 
Magnitude at 
t=6000 seconds 

0° 

3.85 

51.6° 

3.84 

90° 

3.87 
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, there was no significant increase in the attitude errors for a polar orbit. In fact, 
there were not fewer satellites visible when the spacecraft was over a pole than at any 


other position in the orbit, as is shown in Figure 5.11. 



Figure 5.11 Number of Satellites Visible During a Polar Orbit 


5.3.6 Attitude Hold 

Three types of attitude hold were compared: LVLH, inertial, and attitude maneu- 
vers. LVLH maintains the spacecraft orientation with respect to the earth; the x-axis 
points in the direction of velocity, the z-axis points away from the earth, and the y-axis 
is the cross track direction. This attitude corresponds to an aircraft flying nose forward, 
wings level around the earth. When the spacecraft was simulated as being in an inertial 
attitude hold, the initial attitude was the same as the LVLH case, but the inertial atti- 
tude remained constant throughout the run. In the maneuver case, the spacecraft again 
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started in the same initial attitude, and then performed a series of rotations, including 
one about each body axis. These rotations are recorded in Table 5.12. 

The attitude errors were significantly less for the LVLH case than for the inertial 
case, as shown by the summary of results in Table 5.1 1, and in Figure 5.12.The main 

Table 5.11 Performance for Attitude Hold Cases 


Attitude Hold 
Type 

Attitude Error (mrad) 

t=7900 

t=86400 

LVLH 

3.53 

1.96 

Inertial 

5.27 

3.10 

Maneuvers 

2.90 

3.09 



cause of this performance difference is that the rotation about the body y-axis which 
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occurs in the LVLH hold (due to the orbital rotation) allows errors in the baseline coor- 
dinates to be observable, and thus those errors can be estimated.This reasoning is sup- 
ported by the results when maneuvers were performed which included rotations about 
each body axis. The performance results immediately after the conclusion of the 


maneuvers and at the end of the run are included in Table 5.11, and Figure 5.13 shows 



Figure 5.13 Attitude Errors During Maneuvers 


the dramatic decrease in attitude error which occurs during each maneuver. However, 
as Figure 5.12 shows, this decrease is not maintained after the maneuvers have been 
completed. 

When the spacecraft is in an inertial attitude hold, the earth blocks most of the sat- 
ellites during a portion of each orbit. This results in an increase in attitude error, as can 
be seen in Figure 5.14. Figure 5.15 shows the number of satellites visible during this 
time. 
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Number of Satellites Visible Attitude Error Magnitude (mrad) 
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Figure 5.14 Attitude Errors Caused by Loss of Satellite Visibility 



Figure 5.15 Number of Satellites Visible While in an Inertial Attitude Hold 
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Table 5.12 Maneuver History 


Time Period 
(seconds) 

Maneuver Description 

Body Axis of Rotation 

Rotation Rate (deg/hr) 

0- 1000 

None 

0.0 

1000- 1900 

X 

0.1 

1900 - 4000 

None 

0.0 

4000 - 4900 

Y 

0.1 

4900 - 7000 

None 

0.0 

7000 - 7900 

Z 

0.1 

7900 - 86400 

None 

0.0 


5.3.7 Gyroscope Quality 

Two levels of gyroscope quality were modeled, a good quality gyro, which was 
similar to the gyro proposed for the Space Station, and a poor quality gyro. There was 
a marked degradation of performance when the poor quality gyroscopes were used; 
the results are summarized in Table 5.13, and Figure 5.16 shows the attitude errors. 


Table 5.13 Gyroscope Quality Results 


Gyro Quality 

Attitude Error 
Magnitude at 
t = 86400 seconds 
(mrad) 

Good 

1.96 

Poor 

6.89 


The reason that the gyro quality affects the long-term performance of the GPS/inertial 
system is that the baseline coordinates cannot be estimated as well with the bad gyros, 
as shown in Figure 5.17. 
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Gyro Bias (deg/hr) ST Attitude Error Magitude (mrad) 



When the poor quality gyroscopes were used, it was possible to estimate the gyro 


. As Figure 5.18 shows, the estimate of gyro bias was improved by an order of 



Time (s) 


Figure 5.18 Gyro Bias Estimation for a Poor Quality Gyroscope 
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Baseline 1 RMS Errors 


Baseline 2 RMS Errors 




Figure 5.17 Comparison of Baseline Errors for Good and Poor Quality Gyros 
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magnitude over the course of the run, in effect improving a poor quality gyroscope to 
medium quality. The scale factor errors and the misalignments were not significantly 
reduced. 
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Chapter 6 


Conclusions and Suggestions for Future Work 

In this thesis, the performance of a GPS/inertial attitude determination system was 
investigated using a linear covariance analysis. The major error sources of both GPS 
interferometry and gyroscopes were modeled. A new figure of merit (Attitude Dilution 
of Precision, or ADOP) for satellite selection was presented and derived. Finally, the 
results of the linear covariance analysis were presented. 

This chapter discusses the factors which enhance the accuracy of a GPS/inertial 
attitude determination system; the factors that limit accuracy; and finally, the factors 
that are not significant to the level of accuracy achieved. Also, suggestions for future 
research are made. 

6.1 Conclusions 

The accuracy of the attitude determination estimate, given the error sources and 
levels, can be improved in two ways. First, measurements can be taken over a long 
period of time (approximately 15 hours). Second, maneuvers can be performed which 
increase the observability of the errors, and therefore allow for better and faster esti- 
mation. 

However, the estimation of the baseline errors is a limiting factor for the attitude 
accuracy. If the RMS baseline flexures can be reduced, then not only will the attitude 
estimate accuracy improve, but the system will also reach its steady state performance 
more quickly, and lower quality gyros may be employed. 

The simulation demonstrated that several factors which were thought to be poten- 
tially significant were in fact not very significant. For example, among the methods 
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considered, the method used for satellite selection did not significantly affect the 
errors. In addition, there was very little difference in results from the single difference 
vs. double difference measurements. The measurement type chosen will therefore 
depend on other factors, such as computational considerations, the number of satellites 
expected to be visible, and any limitations on the number of filter states which are 
available. 

Finally, it was found that there is no loss of satellite visibility in high inclination 
orbits, and therefore orbital inclination does not significantly affect the attitude errors. 
The amount of multipath correlation also has very little bearing on the size of the 
errors. Increasing the amount of multipath results in only a slight increase in the RMS 
attitude error. 

6.2 Future Work 

In this thesis, an error analysis was conducted using a linear covariance simulation. 
The next step in evaluating the performance of this system is to implement a Monte 
Carlo simulation. In this type of simulation, nonlinear effects could be included in the 
environment model, resulting in a more accurate model. 

The GPS/inertial attitude determination system has the potential to be used for 
“fault detection,” i.e. to reject erroneous GPS measurements, a possibility that was not 
investigated here. Because of the high bandwidth of the gyroscopes, they could easily 
be used to detect cycle slip or an incorrect integer ambiguity. 

In this analysis, baseline errors were modeled as exponentially correlated random 
variables. Because these errors were the most significant, it is important that they be 
properly characterized. In addition, performance may be improved by the develop- 
ment of multipath calibration techniques suitable for space vehicles. 
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Appendix A 


Geometric Dilution of Precision (GDOP) 


In Reference [16], Leick describes GDOR The following is a summary of that dis- 
cussion. 

Generally, the parameters that are important to the positioning problem are 


{x k , y k , z k , dt k } (A.i) 

where x k , y k , and z k are the coordinates of receiver k in an Earth-centered, inertial ref- 
erence system, and dt^ is the clock error of receiver k. The observable most often used 
for positioning is the pseudo range, which is a measure of the distance between the 
receiver and the satellite. The pseudo range equation for receiver k and satellite p is 


Pk = J(x p -x k ) 2 + (y p -y k ) 2 + (z p -z k ) 2 + c • dt k 


= P£ + cdt k 


(A. 2) 


where (x p , y p , z p ) are the coordinates of the satellite, and c is the speed of light. Note 
that Leick’s notation uses a superscript to denote the satellites (“high in the sky”), and 
subscripts to denote a receiver on earth. The design matrix is formed by calculating the 
sensitivity of each pseudo range to each of the parameters listed in Equation (A. 1 ): 
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(A. 4) 


For terrestrial applications, the covariance matrix above is usually transformed to a 
local navigation frame. However, inertial coordinates were the most appropriate for 
the space application considered in this thesis. GDOP is then calculated as root sum of 

T 1 

the diagonal elements of (A A) , so that 


GDOP = -/trace ( A T A) 1 


(A. 5) 



Appendix B 


Derivation of Pitch and Yaw Sensitivity Equations 

The equations presented in the reference [14] for “Elevation Sensitivity” and “Azi- 
muth Sensitivity” of a double difference measurement using satellites i and j are 

Elevation Sensitivity = -cos (ou - a R ) cos^sin^ + sinP^cosf^ (B.l) 
Azimuth Sensitivity = cosp R cosp ij (-cosa ;i sina R + sina^cosotR) (B.2) 
where the variables are: 

a R = antenna baseline azimuth 

P R = antenna baseline elevation 

otj . = azimuth of vector to satellite (s) 

P (j = elevation of vector to satellite (s) 

These equations apply equally to the SD and DD cases, the only difference 
between the two being the vector to the satellite(s). In the SD case it is the LOS vector 
to the satellite, in the DD case, it is the difference of the unit LOS vectors. The satellite 
elevation and azimuth are shown in Figure B.l for both the SD and DD cases. 

In the orbital environment, “azimuth” and “elevation” are not well defined in refer- 
ence to either the baselines or the satellites. Therefore, a local level frame was defined 
with the horizontal plane coincident with the antenna array plane, and the vertical 
direction the same as the antenna “look direction,” which is perpendicular to the 
antenna array plane. 
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Figure B.l Satellite Elevation and Azimuth Angles for DD and SD Cases 

Then, baseline “elevation” was redefined to be baseline “pitch,” or rotation out of 
the antenna array plane. Baseline “azimuth” was redefined to be “yaw,” or rotation in 
the array plane. Since the local level frame is defined by the baseline location and 
direction, ocr and Pr are identically zero. Satellite “azimuth” and “elevation” in the 
local level frame are simply the spherical coordinates of the satellite relative to the 
baseline, as shown in Figure B.l. Therefore, pitch sensitivity can be calculated from 
Equation (B.l) as 

Pitch Sensitivity = -cos (a, - 0) cos^sin (0) + sinPjCos (0) 

= 0 + sinpj ( 1 ) 

= sinB 

1 (B.3) 

and yaw sensitivity can be calculated from Equation (B.2) as 

Yaw Sensitivity = cos (0) cosPj (-cosoesin (0) + sinoccos (0) ) 

= (l)cosp,((0) + sinotj ( 1) ) 

= cos P sin a 

1 (B.4) 
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Appendix C 


The Dynamics and Process Noise Matrices 


Given the state vector x T = 


Vj/ T g T p T y T 8b| SbJ p, p 2 , then the dynamics 
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The continuous process noise matrix is 
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and 


Q P . 



(C.5) 


The discrete noise covariance matrix is 


Qk = 


<ky 
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0 
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0 
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'kbl 


; kb2 


Q k ^ 0 0 

0 Q k ^ 0 
0 0 Q k ^ 


*kp| 


; kp2 


(C.6) 


where, using dt = t k+) - 1^, 


Qk.,< - 


Qarwdt 0 0 | 

0 q arw dt 0 

0 0 qarwdt 


(C.7) 
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and 


Qkpl = Qkp 2 = a p ( 1 - e _2d,/Tp ) 


'kp2 
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Appendix D 


Derivation of the Sensitivity Vectors 

D.l Single Difference Case 

The sensitivity vectors for the baseline errors and path delays can be calculated 
directly from the measurement equation; their derivations follow in Section D. 1 . 1 . The 
sensitivity for the attitude errors is more easily obtained by differencing the actual and 
predicted measurement. The attitude error sensitivity is derived in Section D. 1.2. 

D.1.1 Baseline Errors and Path Delays 

The measurement equation is (from Equation (3.7)) 

SD = X(n + ^) +pd = (uj,) T C^b B + pd 
Since the sensitivity of the baseline errors is the same as the sensitivity to the baseline 
coordinates, hg b is simply the partial derivative of the measurement with respect to b, 
given by 

t>5b = (U|,) T CL ID') 

Likewise, the scalar sensitivity for each path delay is 

h p d = 1 < D - 2 > 

Obviously, a single difference measurement is only sensitive to the baseline errors and 
path delay which are associated with the baseline which was used for the measure- 
ment. 

D.1.2 Attitude Errors 

In general (from Equation (4.13)) 
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^meas b X + Y 


qpred = 

and therefore 

— qmeas q p red 

= -h T 5x + Y 

Omitting measurement noise sources, 

, A<t> 

^Imeas ^ (n+ 2 i )+pd 

= (u p ) T Cg!? B + pd 
and 

q p red = (-Ul,) T CBb B 

Therefore, 

5q = (uJ,) T (C^ B -CBb B )+pd 

By definition 

b B = b B + 5b B 
\|/ B = y B + 8y B 

and 

Cb = (I + fx)C^(I-tj/ B x) 

= C^(I + y B x-y B x) 

= Cjjd + S^x) 


(D.3) 


(D.4) 


(D.5) 


(D.6) 


(D.7) 


(D.8) 


(D.9) 
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where \j/x is the cross product matrix form of the indicated vector. So, 


5q = (Ur) T C g (b B - (I + 6\j/ B x) (b B + 8b B )) +pd 
= (u B ) T Cg (b B - b B - 6b B - 5\|/ B xb B ) + pd 

= (ul,) T C| 1 (-5b B + b B x8ij/ B )+pd (D)0) 

As with the baseline sensitivity, is the same as hg v , so 

bj = -(y^C^dA) (D'l) 

D.2 Double Difference Case 

As stated in Chapter 4, double difference measurement sets are decorrelated before 
being incorporated sequentially. Section D.2.1 describes the sensitivity vector of the 
actual measurement, and Section D.2. 2 describes the transformation which decorre- 
lates the measurements. 

D.2.1 Sensitivity Vectors 

The sensitivity vectors for the baseline and attitude errors are the same as in the 
single difference case, with the exception that the LOS difference, (u^, — u B2 ) , is 
substituted for the LOS, uj*. Therefore, 

fa 6b = (y R 1 Wr2) C b (D.I2) 

and 

bj = — (wki -UR 2 ) Tc B(fa Bx ) (° 13 ) 

The path delays cancel out in double difference measurements, therefore, the path 
delay sensitivity is zero. 
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D.2.2 Measurement Noise Decorrelation 

As explained in Section 4. 1 , R must be diagonal in order to do sequential process- 
ing. Any full square matrix can be transformed to a diagonal one using an eigenvalue/ 
eigenvector decomposition. That is, given any nxn matrix R 


T t RT 



0 


0 


(D. 14) 


where are the eigenvalues of R, and the columns of T are the eigenvectors of R. T is 
called the modal matrix. 

A correlated batch of measurements with the sensitivity vectors forming the col- 
umns of H, and <v y T > = R (not diagonal) is represented by Equation (4.9). An 
uncorrelated set of measurements which contain the same information is given by 


z' = T t z = T t Hx + T t v (D. 15) 

The measurement set defined by Equation (D.15) is uncorrelated because 


(T t v v t T) = T t (v • y T )T 
= T t RT 


(D. 16) 


and T*RT is a diagonal matrix, as shown in Equation (D.14). 

T 

The sensitivity vectors, h , for the measurement z‘ are then the rows of T t H. 
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